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REPORT  SUMMARY 


This  is  the  first  semi-annual  technical  report  describing  the  work 
done  and  the  accomplishments  of  the  research  effort  for  ONR-ARPA  titled, 
"Composite  Design  Synthesis". 

The  research  being  carried  out  deals  with  the  problem  of  design 
synthesis  in  heterogeneous  elasticity.  Design  synthesis  is  defined  as  the 
achievement  of  a  desired  design  criterion,  i.e.,  stress  distribution,  strength- 
to-weight  ratio,  etc.,  by  preselecting  a  stress  or  displacement  pattern  in  a 
stretched  plate  and  then  determining  the  variation  of  the  elastic  moduli  that 
is  required  to  permit  the  desired  effects.  The  accomplishment  of  the  desired 
effects  requires  the  solution  of  the  governing  equations  of  elasticity,  parti¬ 
cularly  the  compatibility  equation,  in  terms  of  preselected  stress  fields  in 
the  body  for  unknown  material  properties  which  are  spatial  functions. 

Since  the  initiation  of  the  program,  progress  has  been  made  in 
several  areas.  The  compatibility  equation  for  generalized  plane  stress  has 
been  generated  in  terms  of  both  rectangular  and  polar  coordinates.  Solutions 
for  material  property  variation  for  the  pressurized  annular  disk  have  been 
generated  wherein  the  stress  condition  imposed  was  one  of  constant  circum¬ 
ferential  stress  through  the  wall  of  the  disk.  This  solution  showed  that 
such  a  stress  distribution  can  be  achieved  if  a  smoothly  varying  modulus  of 
elasticity  is  developed  such  that  it  is  a  given  function  of  the  radius. 
Material  property  relations  of  a  modified  type  of  orthotropy  as  well  as 
isotropy  were  employed  and  curves  plotted  for  a  range  of  material  properties. 

It  was  shown  that  quite  different  distributions  of  the  moduli  are  required  for 
the  externally  pressurized  disk  as  compared  to  the  internally  pressurized  disk. 
These  moduli  functions  were  derived  directly  from  the  compatibility  equation. 
The  work  done  showed  that  it  is  theoretically  possible  to  achieve  a  maximum 
circumferential  stress  in  a  pressurized  disk  that  is  considerably  lower  than 
can  be  achieved  in  the  homogeneous  isotropic  or  orthotropic  case,  and  in  fact, 
for  the  internally  pressurized  disk  this  maximum  stress  can  be  considerably 
lower  than  the  applied  pressure. 

The  governing  equations  of  the  heterogeneous  rotating,  uniform 
thickness,  rotating  annular  disk  have  been  formulated  and  the  solutions  are 
being  determined  numerically.  In  this  case  the  material  density  as  well  as 
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the  \aoduli  can  be  considered  as  spatial  functions,  thus  introducing  another 
degree  of  freedom  that  might  be  employed  by  a  designer. 

A  stress  function  has  been  developed  which  if  implemented  by  the 
correct  material  property  variation  will  eliminate  the  stress  concentration 
found  around  small  holes  in  large  plates.  This  particular  application  has 
great  usefulness  in  the  fabrication  of  structures  that  are  secured  by  rivets 
or  for  structures  which  must  contain  circular  cutouts. 
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SYMBOLS 

Quantity 


a,  b 


c 


A  ,  B  ,C  « C  i  sCa 
o’  o’  o’  1*  2 


E. 

l 

f 

f1(e),f2C0) 
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m 

n 

P 

q 

r 

T 

u,v 

R,X,Y 


inside  and  outside  radius 
material  coefficients 


dimensionless 


shear-modulus  coefficient 


constants  of  integration 


orthctropic  ratio 


modulus  of  elasticity  corresponding  to  the  subscripted 
direction  (i  «=  r,9) 


functions  of  9  only 


modulus  of  rigidity  (  =  G  .  =  - ) 

\  re  a 66J 

acceleration  of  gravity 

radius  ratio  (=b/a) 

orthotropic  ratio  ( =e) 

Poisson's  ratio  in  tangential  direction  (  =v  ) 

U  r 


COS  0 
sin  0 

internal  pressure 
external  pressure 
radius' 

temperature  difference  function 

displacements 

body  forces 


V 


Symbol 


Quantity 


a 

P 

v 

e 

0 

0 

\ 

v 

§ 

P 

0 

Y 

(JO 


coefficient  of  linear  thermal  expansion 
exponent  ( = -^^2) 
material  density 

strain  component  with  one  or  two  subscripts 
body  force  in  tangential  direction 
angular  position 
exponent 

Poisson's  ratio  in  tangential  direction  (  =  v  ) 

u  r 

exponent  (  e  - 

dimensionless  position  ratio  ( = r/b) 
coordinate  system  offset  angle 
stress  function 
rotational  velocity 


INTRODUCTION 


It  has  long  been  recognized  that  structural  elements  composed  of 
composite  materials,  such  as  glass,  boron,  carbon,  or  other  filaments,  embedded 
in  a  suitable  matrix,  such  as  epoxy  or  polyester,  offer  outstanding  strength- 
to-weight  ratios.  The  potential  of  such  materials  is  considered  so  great  that 
material  scientists  and  engineers  believe  that  they  will  form  the  bulk  of  the 
structural  materials  of  the  future.  \Though  the  application  of  such  materials 
has  been  a  growing  part  of  the  stateof  the  art  for  structural  components, 
-articularly  in  the  aerospace  industry,  the  translation  of  the  concept  of 
fibrous  composites  into  a  primary  load  carrying  structure  has  been  and  remains 
a  challenging  process. 

The  present  and  growing  use  of  structural  elements  fabricated  from 
composite  materials  creates  the  need  for  the  development  of  a  rational  analytic 
design  basis,  which  to  a  great  extent  is  presently  non-existent.  This  is  not 
to  be  construed  as  meaning  that  little  or  no  research  on  composite  materials 
has  been  carried  out.  On  the  contrary,  a  large  amount  of  literature  has  been 
generated  dealing  with  both  the  determination  of  the  mechanical  properties  of 
these  materials  and  the  analysis  of  specific  structures  fabricated  from  them. 

In  general,  from  the  microscopic  viewpoint,  research  on  the  mechani¬ 
cal  properties  of  composites  has  dealt  with  the  determination  of  such  properties 
for  materials  having  given  component  elements  ordered  in  fixed  spatial 
relationships.  The  spatial  relation  in  these  cases  might  have  a  high  degree 
of  symmetry,  as  in  long  or  continuous  filament  composites,  or  a  completely 
random  or  homogeneously  disordered  array  as  usually  employed  in  short  carbon 
or  boron  fiber  composites.  Such  research  has  been  directed  towards  the 
creation  of  analytic  or  experimental  methods  of  determining  the  mechanical 
properties  of  composite  structures  in  terms  of  the  known  properties  of  the 
composites'  components,  characteristic  of  this  approach  are  the  works  of 
Sayers  and  Hanley  [1]",  Chen  and  Cheng  [2],  Hill  [3],  and  Gaonkar  [4],  among 
others. 

In  the  analysis  of  structures  composed  of  such  composites,  the 
material  has  generally  been  treated  as  exhibiting  gross,  homogeneous, 
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isotropic  or  anisotropic  mechanical  properties.  These  gross  properties,  when 
entered  into  the  constitutive  equations  defining-  the  material,  have  allowed 
analyses  of  such  structures  through  the  classical  methods  of  the  theory  of 
elasticity,  plates  ar.d  shells,  vibration,  and  others.  Along  these  lines,  the 
concept  of  the  "unidirectional  lamina"  was  introduced  and  utilized  as  the 
"fundamental  unit  of  material"  in  design  and  analysis.  This  procedure 
employs  test  data  obtained  from  a  unidirectional  lamina  as  the  basis  for  the 
design  of  laminated  components  and  structures.  Exemplary  work  done  along 
these  lines  has  been  carried  out  by  Dong  [5],  Tsai  [6,7],  Tsai  and  Azzi  [8], 
Whitney  [9],  and  Whitney  and  Liessa  [10,11],  also  among  others. 

All  of  these  analyses  have  dealt  with  materials  and  structures  that 
have  predetermined  mechanical  and  behavioral  characteristics.  That  is,  once 
the  geometric  array  and  the  constituents  of  the  composite  are  prescribed, 
then  the  mechanical  properties  of  the  material  and  the  response  characteristics 
of  a  structure  fabricated  from  such  a  material  have  been  inherently  established. 
Analysis  will  merely  determine  what  these  properties  and  response  charac¬ 
teristics  are. 

When  dealing  with  composite  materials,  the  analytical  procedures 
discussed  above  appear  to  be  highly  inefficient  in  many  applications.  The 
designer  of  fibrous  composite  structures  is  presented  with  numerous  degrees 
of  freedom  and.  an  opportunity  to  exercise  ingenuity  totally  unavailable  to  him 
with  conventional  materials.  Composites,  whether  filamentary,  fibrous,  or 
sintered  or  fused  metallics,  are  capable  of  being  tailored  to  meet  specific 
requirements.  When  considering  specific  structural  applications  for  such 
materials,  it  would  be  logical  to  assume  that  a  structure  could  be  optimized, 
depending,  of  course,  upon  the  optimization  criterion,  by  varying  the  mechani¬ 
cal  properties  of  the  material  throughout  the  structure.  Further  it  \tfould  be 
logical  to  bypass  analysis  completely  and  define  this  nov7  nonhomogeneous 
structure  by  some  means  of  design  synthesis.  Admittedly,  the  creation  of  a 
design  synthesis  procedure  to  adequately  handle  most  problems  in  structural 
design  is  quite  difficult.  However,  the  concept  of  design  synthesis  to  deter¬ 
mine  the  variation  of  the  mechanical  properties  of  a  material  within  a  structure 
so  as  to  achieve  a  desired  stress  or  deformation  pattern  in  that  structure  is 
one  capable  of  being  developed. 


The  work  described  herein  deals  with  the  first  phase  of  an  effort 
to  develop  a  design  synthesis  methodology  for  composites.  It  is  directed 
specifically  to  the  problem  of  heterogeneous  plane  elasticity. 
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GENERAL  DISCUSSION 


Consider  the  following  question: 

Given  a  plane  elastic  body  with  known  boundary  tractions  and/or 
displacements,  can  the  mechanical  properties  of  the  continuum  be  described 
such  that  an  "arbitrary"  stress  distribution  within  the  body  is  met? 

The  term  "arbitrary"  is  to  be  understood  as  defining  a  family  of 
stress  distributions  that  are  preselected  but  still  conform  to  equilibrium 
requirements  and  boundary  conditions.  To  answer  this  question  we  start  by 
making  the  following  two  basic  assumptions: 

(1)  the  classical  equations  of  linear  elasticity  are  valid  in  this 
application,  and 

(2)  the  mechanical  properties  of  the  continuum  can  be  expressed  as 
spatial  functions. 

Following  from  these  assumptions  the  well  known  governing  relations  of  gener¬ 
alized  plane  stress,  given  in  rectangular  coordinates  are  as  follows: 


Equilibrium  Equations 


da 

_ x 

dx 


da 

XV 

dy 


+  X  <=  0 


!Zz 

dy 


+  Y  f=  0  , 


(1) 


Strain-Displacement  Equations 
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Assuming  that  the  continuum  exhibits  orthotropic  material  proper¬ 
ties  and  neglecting  time  and  strain  rate  effects,,  the  constitutive  equations 
can  be  expressed  as  generalized  Hooke's  Law  as 


where  a  and  a  are  the  coefficients  of  thermal  expansion  and  T  represents 

J.  ^ 

the  temperature  difference  distribution.  The  relations  defined  by  Equations 
(4)  are  valid  when  the  axes  of  the  material  properties  of  the  continuum  are 
coincident  with  the  axes  selected  for  the  differential  equations  of  the 
problem.  If  the  axes  are  not  coincident,  except  for  the  z-axes,  and  the 
other  two  axes  of  the  material  properties  are  rotated  about  the  z-axis  through 
some  angle,  <h,  in  relation  to  the  geometric  axes,  then  more  complicated 
relations  between  the  stresses,  temperature  and  strains  are  developed. 
Lekhnitskii  [ 12 ]  presents  these  relationships  in  some  detail.  For  the  gener¬ 
alized  plane  stress  case  in  point  the  material  coefficients  are  related  to  the 
two  axis  system  by 


where  m  *=  cos  <t,  and  n  =  sin  <  . 

No  generality  will  be  lost  by  continuing  with  the  constitutive 
equation  as  given  by  Equations  (4).  Substituting  Equations  (4)  into  (3), 
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assuming  that  the  a..'s  are  spatially  dependent  and  carrying  out  the 
differentiation  yields: 
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Assume  that  the  body  forces  have  a  potential,  V,  such  that 
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and  choose  Airy's  stress  function,  Y ,  in  the  form 
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Making  the  appropriate  substitutions  into  Equation  (6)  yields 
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and,  of  course,  the  equilibrium  equations  are  met  exactly. 

Notice  that  in  Equations  (4),  (5),  and  (9),  a^2  has  not  been 
equated  to  321*  Normally,  under  the  limit  of  small  displacement  theory 
dealing  with  linearly  elastic  materials  which  are  conservative,  the  material 
coefficient  matrix,  defined  in  Equation  (4),  would  be  symmetric  and  a^2  would 
equal  a^.  However,  some  recent  work  by  Bert  and  Guess  { 13 ] ,  among  others, 
shows  that  there  exists  experimentally  derived  data  which  indicate  that  for 
some  types  of  composite  materials  exhibiting  orthotropic  properties  the 
material  coefficient  matrix  is  not  symmetric  and  a^2  is  not  equal  to  a2^. 

To  limit  the  growing  complexity  of  this  work,  the  material  coefficient  matrix 
will  be  taken  as  symmetric,  at  least  for  the  initial  phase  of  this  effort. 

Equation  (9)  can  also  be  expressed  in  polar  coordinates  as 


follows : 
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with  Y  being  the  stress  function  defined  by 
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Equations  (11)  satisfy  the  equilibrium  equations  formulated  in  polar  coor 
dinates  for  plane  stress  which  are 
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In  order  for  Equations  (11)  to  meet  Equation  (12)  exactly  the  body  functions 

V  and  V  must  be  defined  as 
r  9 


f-  (rV  ) 
9r  v  v7 


R 


99 


0 


(13) 


thus  putting  a  rather  narrow  interpretation  on  the  body  forces. 

In  the  classical  approach  to  the  problem  of  orthotropic  plane 
elasticity,  the  coefficients  a„  are  either  constants,  as  in  the  case  of  the 
homogeneous  condition,  or  as  in  some  rare  cases,  are  special  functions  of 
position.  In  the  first  case  all  the  terms  within  the  second  set  of  large 
braces,  {}•  in  Equations  (9)  and  (10)  are  zero  leaving  the  remainder  of 
these  two  equations  in  the  form  of  the  well-known,  homogeneous,  orthotropic, 
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compatibility  equations  in  terms  of  the  stress  function  Y.  In  the  second 
case  all  the  partial  derivatives  of  the  a..'s  which  appear  in  these  second 
sets  of  braces  are  capable  of  being  evaluated,  resulting  in  extremely  compli¬ 
cated  fourth  order  partial  differential  equations  with  variable  coefficients. 
Proceeding  along  classical  lines,  these  equations  must  be  solved  for  Y  which 
contains  arbitrary  constants  of  integration.  These  constants  are  then  deter¬ 
mined  by  evaluating  Y  in  terms  of  the  stresses  on  the  boundaries. 

Suppose  it  is  assumed  that  the  material  coefficients,  the  a..'s, 
are  unknown  but  that  the  stress  function  Y  is  a  fully  defined  function  of  the 
spatial  coordinates.  That  is,  the  stresses  throughout  the  body  as  well  as  on 
the  boundary  are  known.  In  such  a  case,  Equations  (9)  and  (10)  reduce  to 
second  order  partial  differential  equations  with  variable  coefficients,  in 
terms  of  the  a^'s.  The  solution  of  these  equations  and  the  resultant  deter¬ 
mination  of  the  magnitude  and  distribution  of  the  material  properties  throughout 
the  body  is  defined  as  design  synthesis. 

It  is  believed  that  Equations  (9)  and  (10)  have  not  been  previously 
published.  Bert  [14]  derived  an  equation  similar  to  (10)  wherein  he  reduced 
the  unknown  material  property  coefficients  from  four  to  one.  His  formulation 
is  as  follows: 
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d0 
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J  dr2  Ldr2  ’  rdr  -  r2 
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(14) 
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with  the  thermal  terns  omitted.  Equation  (14)  is  equal  to  Equation  (10) 
with  the  following  identities: 


e  S  all/a22 
e=  a66/2a22 

V  S  "  a12/a22 


(15) 


with  S  being  the  dependent  variable  and  e,  c,  and  v  being  fixed  ratios.  It 

is  anticipated  that  most  real  composite  materials  will  exhibit  material 

property  characteristics  as  defined  by  Equations  (15). 

The  solution  of  equations  such  as  (9)  and  (10)  where  there  are  4 

or  more  independent  a..'s,  even  where  these  a..'s  are  assumed  independent  of 

i  J 

temperature,  is  quite  difficult  but  certainly  not  impossible.  The  simpli¬ 
fying  assumption  made  in  Equations  (15)  leading  to  the  formulation  of  Equation 
(14)  reduces  the  problem  to  only  one  unknown  parameter.  Work  is  progressing 
on  the  application  of  Equations  (9),  (10),  and  (14)  to  the  concept  of  design 
synthesis.  The  following  section  of  this  paper  deals  with  such  applications. 
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APPLICATIONS 


Rotatlonally  Symmetric  Problems 


Example  1:  The  pressurized  annular  disk,  internal  pressure:  Consider  a 
pressurized  annular  disk  as  shown  in  Figure  1.  In  the  case  where  the  ring 
is  isotropic  and  homogeneous,  the  stress  distribution  is  as  developed  by 
Lame  (1852)  and  is  given  by  (c.f.,  Ref.  [15],  p.  60), 


(16) 


where  p  *=  r/b,  and  k  «=  b/a.  These  relations  lead  to  the  following  con¬ 
clusions  : 

(1)  |a0|>|arl  for  all  p  and  all  k 

( 2 )  (n  )  occurs  at  the  inner  boundary  (p  «=  1/k),  and  thus 

v  v  e'max 


^a0  ^max 


>  P 


(17) 


From  the  stresses  so  generated,  which  for  the  homogeneous  ring,  are  quite 
independent  of  the  material  properties,  it  is  clear  that  the  material  is  not 
being  used  efficiently,  particularly  as  the  thickness  ratio,  k,  increases. 

For  the  homogeneous,  ortho tropic  annular  disk  Bienick  et  al.  [16] 
shows  that  the  stresses  are  given  by 


-<£+1>  +  c  *-i 


,  -<f+i)  .  r  f-i 

a0  *=  -c^fr  + 


(18) 


where  f  «=  (a^/a^)1^,  and  C1  and  C2  are  determined  from  the  boundary  con¬ 
ditions,  in  this  case  ar  c  -P  at  r  «=  a  and  ar  c  0  at  r  «=  b.  This  work  and 
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work  done  by  Shaffer  [ i7 1  show  that  as  k  gets  large  the  limit  of  the  stress 
concentration  approaches  the  orthotropy  ratio,  f,-  and  does  so  rapidly.  Both 
Bienick  and  Shaffer  deal  with  the  nonhomogeneous  case,  but  both  assume  a 
special  form  of  nonhomogeneity,  mainly  that 

all  -  all  rX 
a22  *  a22  ^ 
a12  c  a21  *  a12  ^ 

where  \  is  real  and  a^,  a,>,j,  and  a ^  are  constants.  Both  carry  out  optimi¬ 


zation  by  varying  \  and  observing  the  results,  but  no  attempt  was  made  to 
carry  out  design  synthesis  directly. 

Consider  now  that  the  disk  is  composed  of  a  heterogeneous,  ortho¬ 
tropic  material  where  the  material  properties  are  undefined  but  are  functions 
of  the  radius  of  the  disk.  The  equilibrium  equation  for  this  case  reduces  to 


dT+  ; — 0  • 


(19) 


The  compatibility  equation  reduces  to 


fc  (rV 


s  «=  0 


(20) 


and  the  stress  strain  relation  for  an  orthotropic  medium  is 


er  c  allar  +  a12a6 


e9  "  a12CTr  +  a22a0 


(21) 


where  a^,  a^,  and  a£2  are  functions  of  r.  Assuming  a  stress  function 


such  that 
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CT 


r 


CT 


0 


dY 

dr 


(22) 


and  substituting  Equations  (21)  and  (22)  into  Equation  (20)  and  carrying  out 
the  required  differentiation  results  in 


*22*"  +  <a'22+7  a22»'  +  <7  a'l2  *11>*  ’  0 


(23) 


where  the  prime  marks  designate  absolute  derivatives  with  respect  to  r. 

Equation  (23)  is  a  single  differential  equation  with  three  independent 
material  coefficients  as  the  operational  parameters.  When  dealing  with  specific 
materials  they  may  be  found  to  be  completely  independent  or  show  some  type  of 
defined  relationship.  As  a  first  step  we  will  assume  that  there  does  exist  a 
simple  ratio  relationship  between  them  that  is  expressable  as 


CM 

CM 

c3 

*  a22  ’ 

a ' 

22 

“  “'22 

an 

=  \*22> 

a  *  11 

=  kla*  22 

(24) 

al2 

■  k2a22’ 

CM 
!“" 1 

’  k2a'22  ' 

Making  the  appropriate  substitutions  in  Equation  (23)  yields 
k  k 

[y*  +  y  y]  a'22  +  [y  '  '  +  ~  Y '  "  “ 2  y]  a22  *=  0  .  (25) 

r 

Thus  if  Y  is  known  a22  is  fully  defined.  Let  us  suppose  that  for  effective 
material  utilization  it  is  desired  that  cTq  be  constant  throughout  the 
cylinder;  i.e.,  aQ  =  A  .  Integrating  the  second  of  Equations  (22)  gives 

o  O 

Y  «=  A  r  +  B  . 
o  o 

Applying  the  boundary  conditions  that  CTr(a)  **  -  P  and  CTr(b)  =  0  yields 


(26) 
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*  ■  iTTT  [p  -  11 

°r  '  '  FTT  lp  -  11  <27> 


with  p  =  r/b  and  k  «=  b/a  as  before.  It  is  interesting  to  compare  Equations 
(27)  with  Equations  (16)  and  (17).  From  Equations  (27)  it  can  be  seen  that 
if  k  is  greater  than  2,  <7q  becomes  less  than  the  pressure  P,  and  a  ,  which 
equals  P  at  p  «=  l/k,  becomes  the  maximum  normal  stress  in  absolute  value. 
Thus  for  such  a  stress  function  and  geometry  there  is  no  effective  stress 
concentration.  This  is,  of  course,  never  true  for  the  homogeneous  disk. 

With  the  stress  function  now  fully  defined,  Equation  (25)  becomes 


r  k2 

b-,  kib 

[(i+k2)  --f 

"Ja,22+[  r  +  2 

r 

J  a22  c  0  * 

(28) 

which  has  the  solution 


a 


22 


(29) 


where  p  =  r/b 

3  -  -  Vk2 

-  r  Vk2 

5  c  "  L(i-k2)k2 

and  for  the  modified  orthotropic  condition  as  defined  by  Equations  (21)  and 
(24) 


a22  “  a22 


!u 

a22 

fl2 

a22 


1 

E 

E 

E 


e 

0 


r 


v 


0r  * 


(30) 
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The  inverse  of  Equation  (29)  or  Eg(r)  is  shown  plotted  in  Figure  2 
for  various  k^'s  while  k2  was  held  constant  at  0.5.  This  figure  shows  the 
strong  dependency  of  the  modulus  distribution  upon  the  orthotropy  ratio  k^. 

The  ratio  k^  has  a  lesser  effect  as  shown  by  Figure  3.  Here  the  isotropic 
case  (k^  «  1)  is  shown  plotted  for  four  values  of  k2« 

Example  2.  Pressurized  annular  disk,  external  pressure:  Consider  the 
pressurized  annular  disk  as  shown  in  Figure  1  but  with  the  pressure  acting  on 
the  OD  rather  than  the  ID  as  shown.  If  the  stress  criterion  is  retained;  i.e., 
Oq  *=  constant,  then  the  stress  function  becomes 

y  ■=  frr  ti  -  Pkl 


-a-  [i 

k-1  p 


-  k] 


_  ak_ 

k-l 


(31) 


where  q  *=  external  pressure 
k  ■=  b/a 
p  «=  r/b  . 


The  compatibility  equation  (25)  becomes 


k2a*l 

~J  3 


22 


+ 


22 


=  0  . 


(32) 


Equation  (32)  is  the  same  as  Equation  (28)  except  that  a  replaces  b.  The 
solution  of  Equation  (32)  is 


(33) 


where  0,  §,  and  p  are  defined  in  Equation  (29).  Comparing  Equations  (33) 
and  (29)  we  see  that  the>  are  not  the  same,  differing  by  the  quantity  1/k 


FIGURE  2.  MATERIAL  PROPERTY  VARIATION  FOR  PRESSURIZED  DISK 
WITH  oa  “  CONSTANT,  INTERNAL  PRESSURE 


20 


within  the  second  factor.  Thus  the  distribution  of  the  material  parameters 
through  the  thickness  of  the  disk  must  be  different  as  compared  to  the 
internally  pressurized  disk.  This  is  shown  in  Figures  4  and  5.  Also  notice 
that  for  the  internally  pressurized  disk  the  hoop  stress  cr  tends  to  zero  as 

U 

k  tends  to  infinity  while  the  limit  on  a«  for  the  externally  pressurized  disk 

U 

is  q,  the  pressure.  Thus  for  the  externally  pressurized  case  the  material  is 
not  being  as  effectively  utilized  and  perhaps  some  other  criterion  might  more 
suitably  apply. 

If  the  annulus  has  both  internal  and  extern ~ .  pressure  and  the  same 
stress  criterion  is  applied,  i.e. ,  Oq  -  constant,  then 


Pb 

k-1 

[p-i] 

-  fe  t1  - pk] 

P 

r1  -  ii 

.  _s_  [ 1  .  kl 

"k-1 

Lp  J 

k-i  Lp 

1  ! 
k-1 

P  -  qk] 

» 

and  the  material  property  variation  is  given  by 


22 


[?♦ 


(34) 


(35) 


Example  3:  The  rotating  disk:  Consider  the  rotating,  uniform  thickness 
annular  disk  as  shown  in  Figure  6.  The  equilibrium  equation  governing  this 
case  is 

d^  (rCTP  "  CT0  +  g  =  °  '  (36) 


The  compatibility  equation  in  terms  of  strain  and  the  stress-strain  relations 
for  an  orthotropic  material  are  given  by  Equations  (20)  and  (21),  respectively. 
Substituting  Equation  (21)  into  (20)  results  in 


12  r 


12' 


■b  3  22ct0  a22?  q 


12  “11' 


22  12'  0 


tej* 


(37) 


FIGURE  4.  MATERIAL  PROPERTY 
FOR  rr  *=  CONSTANT 
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i 


rtCURK  2*  MATERIAL  FtUlFERTT'  VARJATTOH  FOR  FTlESSDRITED  ANNULAR  OISK 
FOT  =  CONSTANT,  EKT&H NA L  PRESSURE  (ISOTROPIC  CASE) 
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Assume  the  stress  function 


Y '  +  —  tu2r2 
g 


(38) 


It  is  possible  that  the  material  density,  y>  could  be  a  function  of  the 
radius.  At  this  stage  such  a  condition  introduces  an  unnecessary  complica¬ 
tion  and  will  not  be  considered.  This  being  the  case  the  governing 
compatibility  equation  becomes 


T"(a22)  +T’(a’22  +  £  a22>  +  * 


12 


,  V  2  2 
+  -*•  u)  r 
g 


r)  =  0  • 


(39) 


Again,  for  simplicity,  the  modified  orthotropy  relations  as  defined  by 
Equations  (24)  will  be  adapted  making  Equation  (39)  take  the  form 


,22[’" 


,2  w  ,y  2  2 

+  —  Y  +—  m  r 
r  gg 


;]+a22[T" 


Y' 


y  2 

— m~  +  —  to  r 
r2  8 


(3-k2)] 


(40) 


The  stress  criterion  will  be  the  same  as  before,  i.e.,  cTq  c 
constant.  Using  this  criterion  and  operating  on  the  second  of  Equations  (38) 
together  with  the  assumed  boundary  conditions  that  ar(a)  =  ar(b)  =  0  yields 


(41) 


25 


where  p  *=  r/b 
k  *  b/e 

v  c=  cub  <=  tip  velocity 

U)  *  rotational  velocity,  radians 

g  *  acceleration  of  gravity. 

Here  we  note  that  in  the  limits,  when  k  -*  1 

Vv2 

ae  m~ 


(42) 


which  is  the  stress  in  a  rotating  thin  ring,  and  when  k  -»  00  (i.e.,  when 
very  small) 


Yv 

re  *  3T  ’ 


(43) 


which  is  smaller  than  exists  for  the  isotropic,  homogeneous  case  by  the 
ratio 


CT0 

(— ) 

'°0O  max 


9+v 


(44) 


where  CT  *=  hoop  stress  for  heterogeneous  case 
01 

a  «=  hoop  stress  for  isotropic,  homogeneous  case 
9o 

v  *=  Poisson's  ratio  for  isotropic,  homogeneous  case  . 


Substituting  Equation  (41)  into  Equation  (40)  and  carrying  out  the 
required  differentiation  yields 


22 


PAr;  +  br  _±_D  1  „  0 

Fr4  +  Gr2  +  HrJ  22 


(45) 


where 


A 


(kL  -  3k2) 


r,3  3-. 

<*  -  ki>  [-rrr] 
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D  •=  k^(ab)  (b  +  a) 
F  «=  k2 


G 


(1.  +  k2) 


rk3  3-, 

fb  -  a 

L  b  -  a  J 


H  *=  k2(ab)  (b  +  a) 


Equation  (45)  is  not  readily  solvable  in  closed  form.  Solutions  are  presently 
being  developed  utilizing  a  finite  difference  approximation  technique. 


Non-Symmetric  Problems 


Example  4:  Small  hole  in  an  infinite  plate.  Figure  7  represents  a  small 
hole  in  an  infinite  plate  which  is  subjected  to  a  uniform  tensile  stress,  P, 
in  the  x-direction.  For  the  homogeneous,  isotropic  condition,  the  stress 
distribution  around  the  hole  is  well  known  as  given  by  Timoshenko  [15]  as 


The  maximum  stress  occurs  at  r  «=  a,  9  *=  (rr/2, 3tt/2)  ,  and  is 


max 


(aq)  , 

9  r*=0,9=n/2 


3P 


For  the  homogeneous,  anisotropic  condition  work  by  Green  and  Zerna  [18], 
Hearmon  [ 19 ] ,  Savin  [20],  Leckhniskii  [12],  and  (as  directly  applied  to 
composites)  by  Greszczuk  [2l],  shows  that  the  maximum  stress  at  the  hole  is 
always  greater  than  for  the  isotropic  case  and  can  reach  values  as  high  as  9P. 


FIGURE  7.  SMALL  HOLE  IN  INFINITE  PLATE  SUBJECTED  TO 
UNIFORM,  UNIAXIAL  TENSILE  STRESS  FIELD 
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Referring  to  Figure  7,  consider  the  portion  of  the  plate  within 
a  concentric  circle  of  radius  b,  large  in  comparison  with  a.  It  can  safely 
be  assumed  that  the  stresses  at  radius  b  are  effectively  the  same  as  in  a 
plate  without  the  hole  and  can  be  given  by 


(<yj  (1  +  cos  2  9) 

r  r=b  2 

,  (47) 

<<Wri=b~-  2  P  Sin26  • 


From  Equation  (46)  it  can  be  seen  that 

(<Tft)  «  7  P  (1  -  cos  20), 
°  r*=b  ^ 


It  seems  reasonable,  then,  to  choose  a  stress  criterion  for  the  plate 


CTe 


function  of  6 


cos  2  6) , 


(48) 


From  the  second  of  Equations  (11),  with  V  equal  to  zero  the  stress  function 
becomes 

Y  ■  \  Pr2  (1  -  cos  2  9)  +  f][(P)r  +  f2(0)  .  (49) 


where  f1(9)  and  f20)  are  functions  of  0  only.  Applying  the  first  and  third 
of  Equations  (11)  to  Y  results  in 


1  SY  ,  1  52Y 
r  9r  r2  S92  =  ^ 


-  ar9  "2  P  Sin20  "  h 

r 


df2(0) 
d0  “ 


d2f  (9)  .  d2f  (0) 

|  P  (1  +  cos  2  0)  +  i  2~+7 

rd0  r  d0 


(50) 

fL(e) 


Applying  the  boundary  conditions 


(ar) 


r*=a 


(crr9) 


r>=a 


0 


yields 


(51) 
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f2(0)  *  “  ~4~  cos  2  0 

d2f,  (9)  p  r 

- +  f t (9)  <=  -  |  [1  +  3  cos  20]  . 

d92  1 


Choosing  as  a  particular  solution  to  the  second  of  Equations  (51) 

f1(9)  “  C1  +  C2  cos  2  0 

yields 

f.,  (0)  *=  -  [l  -  cos20] 

JL  Z 

and  results  in 


Y 

a 


■|  j^(r2  -  2ar)  -  (r2  -  a2)  cos  2  0 j- 


K1-?) 


1  -  cos  20  +  2 


(l  -^)cos  2@] 


a©  ■=  2  (1  "  cos  2  0) 


CTr0  " 


-IM 


sin  2  0 


(52) 


(53) 


The  stress  function  Y  as  defined  by  the  first  of  Equations  (53)  has  been 
applied  to  Equation  (10)  with  the  body  functions  and  the  tetnperature  difference 
distribution  taken  as  zero.  Further,  Equations  (24)  have  been  utilized  to 
initially  simplify  the  problem  and  numerical  solutions  are  being  developed. 
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SUMMARY 


Mathematical  design  synthesis  has  been  shorn  to  be  possible  in 
certain  specific  applications.  The  selection  of  a  design  criterion,  in  the 
cases  discussed,  one  dealing  with  stress  distribution,  and  the  development 
of  the  material  distribution  within  a  plane  body  such  that  compatibility  is 
satisfied,  appears  to  be  a  rational  basis  of  design  for  composite  materials. 
It  is  to  be  shown  that  this  concept  of  design  synthesis  can  be  applied  to 
the  fabrication  of  real  composite  structures. 
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